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Topological insulators 



■ llh are characterized by an insulating bulk with a finite band gap and 
conducting edge or surface states, where charge carriers are protected against backscattering. These 
states give rise to the quantum spin Hall effect |2] without an external magnetic field, where electrons 
with opposite spins have opposite momentum at a given edge. The surface energy spectrum of a three- 
dimensional topological insulator yl, la] is made up by an odd number of Dirac cones with the spin 
locked to the momentum. The long-sought yet elusive Majorana fermion is predicted to arise 



from a combination of a superconductor and a topological insulator 11131 - 11511 . An essential step in the 
hunt for this emergent particle is the unequivocal observation of supercurrent in a topological phase. 
Here, we present the first measurement of a Josephson supercurrent through a topological insulator. 
Direct evidence for Josephson supercurrents in superconductor (Nb) - topological insulator (Bi2Te3) 
- superconductor e-beam fabricated junctions is provided by the observation of clear Shapiro steps 
under microwave irradiation, and a Fraunhofer-type dependence of the critical current on magnetic 
field. The dependence of the critical current on temperature and length shows that the junctions 
are in the ballistic limit. Shubnikov-de Haas oscillations in magnetic fields up to 30 T reveal a topo- 
logically non-trivial two-dimensional surface state. We argue that the ballistic Josephson current is 
hosted by this surface state despite the fact that the normal state transport is dominated by diffusive 
bulk conductivity. The lateral Nb-Bi2Te3-Nb junctions hence provide prospects for the realization of 
devices supporting Majorana fermions 
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Almost simultaneous to the theoretical prediction of 3D topological insulator states in the bis- 
muth compounds Bii-^^Sba;, Bi2Se3 and Bi2Te3 [1], angle-resolved photoelectron spectroscopy 



indeed revea 



compounds [|4|, [Z, 



ed a 



inear dispersion and a helical structure of the Dirac cone at the surface of these 



, ISQ. Soon after, the topological nature of the surface states was confirmed by 
transport studies such as Aharonov-Bohm oscillations in Bi2Se3 nanoribbons Shubnikov-de- 
Haas oscillations in Bi2Te3 lll7H2Jj], scanning tunneling spectroscopy of the square-root magnetic 



field dependence of the Landau level spacing jllll and interference effects resulting from impurity 
scattering jlo|. Now the existence of the topological surface states has been established it is time 
to study the interaction with other materials. Efforts have been made to contact a topological insu- 
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lator to a superconductor in the search for the Majorana fermion. First attempts 122 
the presence of a supercurrent through a topological insulator, although no Josephson effects have 
been observed. Li this Letter we demonstrate that it is possible to induce a Josephson supercur- 
rent in the surface state of the topological insulator Bi2Te3, as evidenced by the two hallmarks of 
the Josephson effect: a Fraunhofer-like magnetic field modulation of the critical current and the 




FIG. 1: E-beam lithographically defined Nb electrodes on exfoliated Bi2Te3. (a) Atomic force mi- 
croscopy image of an exfoliated Bi2Te3 surface. The step edges are 1.0 nm high, corresponding to the 
Bi2Te3 quintuple unit cell. These nanometer flat surfaces span an area up to 50 x 50 jxm^. (b) Scanning 
electron microscopy image of Nb-Bi2Te3-Nb Josephson junctions. The Nb superconductor strips are de- 
fined by e-beam lithography on a 200 nm thick exfoliated Bi2Te3 flake and connected to large Nb electrodes 
fabricated by photolithography. The junctions have a width of 500 nm and a length of (from left to right) 50, 
100, 150, 200, 250 and 300 nm. A supercurrent has been identified at 1.6 K in all junctions up to 250 nm. 
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appearance of Shapiro steps upon microwave irradiation, resulting from the dc and ac Josephson 
effects, respectively. 

We have fabricated poly crystalline Bi^Te-^ samples with a common c-axis orientation using the 



Czochralski method as described elsewhere 1241] . From these, flakes were produced by mechanical 
cleaving, proven to be a powerful tool for the fabrication of a few quintuple layer thin Bi2Te3 flakes 



n25\\ . Figure[T^ shows an atomic force microscopy image of the surface of typical exfoliated Bi2Te3 
flakes. We have performed experiments on large flakes with atomically flat terraces, extending 
over several micrometers. The Bi2Te3 quintuple layer step edges are clearly visible as shown in 
Fig. [I^. The exfoliated flakes have a thickness of 20 nm — 2 \xm. Van der Waals forces bind 
the flakes to the Si-substrate, enabling subsequent sputtering of Nb electrodes by standard lift-off 
photolithography. E-beam lithography is used in a second step to define smaller Nb nanostructures 
by lift-off sputter deposition on top of the flake, connected to the larger Nb electrodes fabricated 
with photolithography. A scanning electron microscopy image of such a device is given in Fig. 

To characterize the electronic properties of our devices and to track them down to the sur- 
face states of a topological insulator we have performed magnetotransport experiments on Bi2Te3 
flakes in a Van der Pauw geometry in fields up to 30 T. From the low field Hall coefficient 
Rh = 7.5 X 10^^ fim/T and the zero-field resistivity p = 300 |ji7cm we deduce a bulk conductivity 
channel with an n-type carrier concentration n = 8.3 x 10^^ cm^'^ and mobility /i=250 cm^A'^s 
corresponding to an electron mean free path /g = 22 nm for a parabolic band. This low mobility 
would not allow the observation of quantum oscillations since the necessary condition, pB ^ 1, 
is far from being fulfilled in the magnetic fields applied. As shown in Fig. |2k, we nevertheless 
observe clear Shubnikov-de Haas oscillations oscillations in the resistance implying another con- 
duction channel with larger mobility. From the angle dependence of the position of the peaks we 
conclude that this additional channel is of two-dimensional nature. 

The oscillatory contribution to the resistance of such a two-dimensional system can be written 



asm 



X . f2'KEF 
Rxx oc — rre" cos — + n + ^B] (1) 



sinh A \ hcOc 

where, Uc is the cyclotron frequency, Ep is the Fermi energy, A = 27<-'^ kBT/hcUc, \d = 
2Tc'^kBTj:,/hujc with the Dingle temperature, and cps is the Berry phase. When plotting the 
1 / i?-positions of the minima and maxima of R^x (Fig. [2]1) as a function of the Landau level index 
it already becomes clear that the corresponding 1 / i?-positions do not extrapolate to = as one 



would expect from an ordinary system, but is rather shifted by ^ implying that = tt, consistent 
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with a half-filled zeroth Landau level that is present in a topological surface state 1117 - 
similar to graphene [28, From the slope in Fig. [2li we infer the carrier concentration of this 
surface state oin = 1.2 x 10^^ cm^^; an effective mass m* = 0.16 mo and a Dingle temperature 
of Td = 1.65 K (corresponding to a mobility ji = 8300 cm^/Vs) follow from the temperature and 
field dependence of the Shubnikov-de Haas oscillations amplitudes (Figs.|2]3 and [2]:). This results 
insLVp = 1.4 X 10^ m/s with /g = 105 nm. 

In high magnetic fields two beating frequencies are present, as Xiong et al. n2(M also observed 
in their study on the topological surface states of Bi^Te^Se. A double frequency has also been 
observed in the 3D topological insulator strained HgTe pTI] . where it was concluded to originate 
from the top and bottom topological surface states, with different carrier densities due to the dif- 
ference in electrostatic environment (substrate and vacuum). As can be seen in Fig. |2^, the peak 
splitting becomes more pronounced with increasing parallel field, though it is already visible in 
the derivative of the magnetoresistance at smaller angles. 

The superconducting Nb-Bi2Te3-Nb junctions, as depicted in Fig.IIh, show Shubnikov-de Haas 
oscillations at high magnetic fields with the same frequency, revealing the presence of topological 
surface states in these devices as well. The critical current of the Nb electrodes on top of the Bi2Te3 
is larger than 30 mA for a strip of 500 nm width and 70 nm height. The junctions have a metallic 
temperature dependence, but upon cooling below 6.5 K, the resistance vanishes completely and a 
supercurrent is observed. Figure[3^ shows a current- voltage characteristic at 1 .6 K showing a clear 
supercurrent of Jc = 18 fiA. To test whether this is a true Josephson supercurrent we performed 
phase- sensitive experiments. Figure [3]1 shows the Fraunhofer pattern from the critical current 
dependence on the magnetic field due to the dc Josephson effect. The smaller sidelobes follow 
from the small ratio between the length and the width of the junction, as predicted by OOll . Due 
to the ac Josephson effect, microwave irradiation should lead to Shapiro steps. Figure [3^ shows 
the JV^-curve upon 10.0 GHz microwave irradiation. Clear steps are observed at integer values of 
V = h/2ef-Rp = 20.7 |aV. A colormap of the conductance dependence on the power and the current 
is shown in Fig. [3h, visualizing the evolution of the Shapiro steps. In Fig. [3]: we have plotted the 
power dependence of the thirst three steps, following the expected Bessel function dependence. 

Now that we have confirmed the Josephson nature of the devices, we can discuss in which 
band of the Bi2Te3 the proximity effect is induced. Figure |4] shows the temperature and length 
dependence of the critical current of the junctions. While the scaling with length can be described 
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by diffusive as well as ballistic transport theory (see Suppl. Information for details), the temper- 
ature dependence of the junctions is clearly far from the diffusive limit and can only be fitted by 
the Eilenberger theory for ballistic junctions 113 iL 13211 . The bulk mean free path, /g = 22 nm, is 
too small to explain the ballistic nature of the supercurrent in the junctions. However, the surface 
conduction band has a larger mean free path = 105 nm. Noting that the mean free path obtained 



from Shubnikov-de Haas oscillations is usually an underestimate [13311 . we conclude that the bal- 
listic Josephson supercurrent is carried by the topological surface states of Bi2Te3. As we discuss 
in more detail in the Suppl. Information, the small IcRn product (reduced to about 2%) follows 
automaticaly from the large bulk shunt. The clean limit fit with Tc = 6.5 K gives a coherence 
length ^ = 2nkBT ~ ^-^ imply iiig th^t = 1-0 x 10^ m/s. This Fermi velocity 

is comparable to the value obtained from the magnetoresistance oscillations of the surface states 
(1.4 X 10^ m/s). 

The possibility of contacting many electrodes on the exfoliated Bi2Te3 flakes combined with 
the reproducibility as observed in the length dependence of the junctions creates an interesting 
basis for experimental investigations of the still increasing number of theoretical proposals on 
topological insulator - superconductor devices il3l- ll5|l . Interestingly, whereas bulk transport tends 
to obscure the observation of surface state transport in the normal state, the supercurrent is found 
to be carried mainly by the surface states. The realization of supercurrents through topological 
surface states is an important step towards the detection of Majorana fermions. 

Methods 

By mechanical exfoliation on polycrystalline Bi2Te3 samples, flakes of size up to 50 x 50 |jm^ 
and thicknesses of 20 nm up to 2 |am were transferred to Si substrates. Sputter deposited 
Nb(200nm)/Pd(5nm) electrodes where defined by lift off techniques using optical photolithog- 
raphy with image reversal photoresist. The Pd layer prevents the Nb from oxidation and the 
thickness of the Nb being on the order of the flake thickness ensures a superconducting contact 
(Ic > 30 mA) between the Nb on the substrate and the flake. The Nb(70 nm) junctions where 
defined with lift off using e-beam and have a superconducting contact {Ic > 30 mA) with the Nb 
electrodes. To obtain transparent contacts, the Bi2Te3 surface is Ar sputter etched in situ prior to 
Nb deposition, while the Bi2Te3 in between the Nb forming the actual junction is unaffected due to 
the lift off technique. The supercurrent is investigated in shielded cryostats with dedicated home 
build equipment and the angle dependent magnetoresistance is measured at the High Field Magnet 
Laboratory (Nijmegen) in a 33 T magnet. 
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- SUPPLEMENTARY INFORMATION - 



Theoretical models for Josephson current 

Any hybrid structure containing superconductors can be described on the basis of the Gor'kov 
equations Im]- In practice, these equations are typically simplified by a quasi-classical approxi- 
mation, which is justified as long as the Fermi- wavelength is much smaller than other length scales 
in the problem. For superconductor - normal metal - superconductor (SNS) Josephson junctions 
Eilenberger quasi-classical equations [|350 are used when the elastic mean free path /g is larger 
than the length L and the coherence length ^. The electronic transport in this clean limit is ballistic 
across the N layer. In the dirty limit of /g ^ L, ^, transport is diffusive and the Usadel equations 
ll36ll are used. When the transparency between the S and N layers is not unity, additional insulating 
barriers (I) are typically included. 

Eilenberger theory fit 

The clean limit theory on the basis of the Gor'ko v eq uations for short SINIS junctions with 
arbitrary barrier transparency D [|3lll was generalized [1321 . 13711 for arbitrar y ju nction length on the 
basis of Eilenberger equations. The supercurrent density J is found to be [32] 

tl(/i)t2(/i) 



J = — efcpfcsTsinx / fidjj,-^ 



where /i = K/kp, ti,2 = -Di,2/(2 - £'1,2), and 

2uJr,L 



Q 



titacosx + 1 + {tit2 + 1) 



A2 



cosh 



+ [ti + t2) . „ smh 



(2) 



A2 



fihvp 



(i-tl) {i-tl) 
(3) 



where the Matsubara frequency is given by a;„ = 2TrkBT{2n + 1), and = ^Jbj\^ /S?. A is the 
gap in the S electrodes, x th^ phase difference across the junction, while vp is the Fermi velocity 
of the normal metal interlayer. The integral runs over all trajectory directions and can be adjusted 
to actual junction geometries. 
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Eq. ^ was evaluated as function of junction length and fitted to the measured critical current 
density. Since the prefactors in Eq. ^ implicitly contain the normal state resistance, which is 
not known for our junctions due to the bulk shunt, we left the overall scale of J free in the fit. 
Subsequently, the best fit to the data at 1 .6 K was obtained for ^ = = 75 nm. It was found 

numerically that the value of the barrier transparencies in the symmetric case had no influence on 
the fitting value for ^. 

The temperature dependence of the critical current was calculated using Eq. Q and the ob- 
tained coherence length. The fit to the measured data is excellent, considering that only the overall 
scale of J was free in this case. 

In fact, the overall scaling factor of the critical current in Eq. (1) can be estimated as well. 
The transparency of the interfaces between the topological insulator and the superconductor are 
important in this respect. The high transparency of our interfaces can be determined from the I{V) 
characteristic. The excess current in the I{V) characteristic is about 67% of the critical current Ic- 



In the Blonder- Tinkham-Klapwijk model ||40D, this gives a barrier strength of about Z = 0.6. For 
these high transparencies, at the lowest temperatures and for the 50 nm junction, Eq. (1) provides 
an IcRn product of the order of 1-2 mV. However the 3D bulk shunt will strongly reduce this 
value by decreasing _R„. From the Shubnikov de Haas oscillations we can estimate the surface 
to bulk resistance ratio. The surface resistance can be found through the amplitude of Rsdh = 
Rr e~^°. With -r4-Te^^^ = 0.77 we estimate a surface resistance of about 2% of the total 

^ Sinn A Sinn A ' 

resistance. Thus, the estimated surface resistance is approximately 29 fi, which together with Ic = 
32 |jA results in IcRn = 1 mV, agreeing with the Eilenberger model. This quantitative agreement 
between model and measurements underlines the conclusion that supercurrent is flowing through 
the ballistic channels of the topological surfaces states, shunted by a normal state bulk conduction. 
Usadel theory fit 



The Usadel equation ll36ll for the S and N layers in a diffusive SNS junction can be written as 



where $ is defined in terms of the normal Green's function G and the anomalous Green's function 
F by $G = cOnF. The normalization condition FF* + = 1 then gives 

Gs,N = I (5) 
IujI + ^s,n%n 



The coherence length is given by ,^ = \i t where D = vple/^ is the diffusion constant. 



The pair potentials As,n are given by 

As,N^n- h 2TckBT > = (6) 

' w„>0 

In the dirty limit, Zaitsev's effective boundary conditions for quasi-classical Green's functions 



were simplified by Kupriyanov and Lukichev 



39|] . When $ and G are found using these boundary 



conditions, finally the supercurrent density can be obtained from 

J = ^Im V -f^N-r^N, (7) 

where Pn is the N layer resistivity. 

For junctions with arbitrary length and arbitrary barrier transparency, no analytical expressions 
exist for the Green's functions. Therefore a numerical code was used to fit the data. In the effective 
boundary conditions [39], two parameters play a role, 7 = ^^^^ and 7^ = ^ {^'jf')' where the 
average of the transparencies takes place over all trajectory angles. For Nb as S electrode and 
Bi2Te3 as N interlayer, 7 <C 1 because of the lower resistivity p of Nb as compared to Bi2Te3. 
The junctions transparency is not known a priori, but from the voltage drop over a barrier in the 
normal state, as well as the large amount of excess current (more than 50% of the critical current) 
in the current-voltage characteristics of the superconducting state, a conservative estimate gives 
D > 0.5, which implies 7^ < 1. Within this parameter range (or even outside the range) no 
consistent fit could be made to the data. Figure 6a shows the fit to the temperature dependence 



of the critical current for 7 = 0.1, 7^ = 1 and ^(Tc) = y 2Jk^T =21 nm, the latter value as 
obtained from fitting the length dependence of the junction. 



* Also at Leiden Institute of Physics, Leiden University, P.O. Box 9506, 2300 RA Leiden, The Nether- 
lands 

[1] Bemevig, B.A. & Hughes, T.L. & Zhang, S.C. Quantum spin Hall effect and topological phase transi- 
tion in HgTe quantum wells. Science 314, 1757-1761 (2006). 

[2] Konig, M. et al. Quantum spin Hall insulator state in HgTe quantum wells. Science 318, 766-170 
(2007). 

[3] Fu, L. & Kane, C.L. & Mele, E.J. Topological insulators in three dimensions. Phys. Rev. Lett. 98, 
106803 (2007). 



[4] Hsieh, D. et al. A topological Dirac insulator in a quantum spin Hall phase. Nature 452, 970-974 
(2008). 

[5] Zhang, H. et al. Topological insulators Bi2Se3, Bi2Te3 and Sb2Te3 with a single Dirac cone on the 
surface. Nature Phys. 5 438-442 (2009). 

[6] Qi, X.L. & Rundong, L. & Zang, J. & Zhang, S.C. Inducing a magnetic monopole with topological 
surface states. Science 323, 1184-1187 (2009). 

[7] Chen, Y.L. et al. Experimental reaUzation of a three-dimensional topological insulator, Bi2Te3. Sci- 
ence 325, 178-181 (2009). 

[8] Hsieh, D. et al. A tunable topological insulator in the spin helical Dirac transport regime. Nature 460, 
1101-1105 (2008). 

[9] Peng, H. et al. Aharonov-Bohm interference in topological insulator nanoribbons. Nature Mat. 9, 
225-229 (2010). 

[10] Zhang, T. et al. Experimental demonstration of topological surface states protected by time-reversal 

symmetry. Phys. Rev. Lett. 103, 266803 (2009). 
[11] Cheng, P. et al. Landau quantization of topological surface states in Bi2Se3. Phys. Rev. Lett. 105 

076801 (2010). 

[12] Majorana, E. Teoria simmetrica dellelettrone e del positrone. Nuovo Cimento 14, 171-184 (1937). 
[13] Fu, L. & Kane, C.L & Mele, E.J. Superconducting proximity effect and Majorana fermions at the 

surface of a topological insulator. Phys. Rev. Lett. 100, 096407 (2008). 
[14] Nilsson, J. & Akhmerov, A.R. & Beenakker, C.W.J. Splitting of a Cooper pair by a pair of Majorana 

bound states. Phys. Rev. Lett. 101, 120403 (2008). 
[15] Tanaka, Y. & Yokoyama, T. & Nagaosa, N. Manipulation of the Majorana fermion, Andreev reflection, 

and Josephson current on topological insulators. Phys. Rev. Lett. 103, 107002 (2009). 
[16] Hasan, M.Z. & Kane, C.L. Colloquium: topological insulators. Rev. Mod. Phys, 82, 3045-3067 (2010). 
[17] Qu, D.X. & Hor, YS. & Xiong, J. & Cava, R.J. & Ong, N.P. Quantum oscillations and Hall anomaly 

of surface states in the topological insulator Bi2Te3. Science 329, 821-824 (2010). 
[18] Xiu, F. etal. Manipulating surface states in topological insulator nanoribbons. Nat. Nano. 6, 216-221 

(2011). 

[19] Analytis, J.G. et al. Two-dimensional surface state in the quantum limit of a topological insulator. Nat. 
Phys. 6 960-964 (2010). 

[20] Xiong, J. & Petersen, A.C. & Qu, D. & Cava, R.J. & Ong, N.P. Quantum oscillations in a topological 



insulator Bi2Te2Se2 with large bulk resistivity (eQcmV E^ivillOl.lBfSV l (2011). 
[21] Taskin, A.A. & Ren, Z. & Sasaki, S. & Segawa, K. & Ando, Y. Observation of Dirac holes and 

electrons in a topological insulator. Phys. Rev. Lett 107, 016801 (2011). 
[22] Zhang, D. et al. Observation of the superconducting proximity effect and possible evidence for Pearl 

vortices in a candidate topological insulator. Phys. Rev. B 84, 165120 (2011). 
[23] Sacepe, B. et al. Gate-tuned norrmal and superconducting transport at the surface of a topological 

insulator. Nature Comm. 2, 575 (201 1). 
[24] Li, A.H. et al. Electronic structure and thermoelectric properties of Bi2Te3 crystals and graphene- 

doped Bi2Te3. Thin Sol. Films. 518 57-60 (2010). 
[25] Teweldebrhan, D. & Goyal, V. & Rahman, M. & Balandin, A.A. Atomically-thin crystalline films and 

ribbons of bismuth telluride. Appl. Phys. Lett. 96, 053107 (2010). 
[26] Lifshitz, I.M. & Kosevich, A.M. Theory of magnetic susceptibility in metals at low temperatures. Sov. 

Phys. JETP2, 636-645 (1956). 
[27] Briine, C. et al, Quantum Hall effect from the topological surface states of strained bulk HgTe. Phys. 

Rev. Lett. 106 126803 (2011). 
[28] Novoselov, K.S. et al. Two-dimensional gas of massless Dirac fermions in graphene. Nature 438, 

197-200 (2005). 

[29] Zhang, Y. & Tan, YW. & Stormer, H.L. & Kim, R Experimental observation of the quantum Hall 

effect and Berry's phase in graphene. Nature 438, 201-204 (2005). 
[30] Barzykin, V. & Zagoskin, A.M. Coherent transport and nonlocality in mesoscopic SNS junctions: 

anomalous magnetic interference patterns. Superlat. and Microstruct. 25, 797-807 (1999). 
[31] Brinkman, A. & Golubov, A.A. Coherence effects in double-barrier Josephson junctions. Phys. Rev. B 

61, 11297 (2000). 

[32] Galaktionov, A.V. & Zaikin, A.D. Quantum interference and supercurrent in multiple-barrier proxim- 
ity structures. Phys. Rev B 65, 184507 (2002). 

[33] Dietl, T. Dingle temperature in HgSe. /. Phys. CoUoques 39, 1081-1083 (1978). 

[34] Gor'kov, L.R On the energy spectrum of superconductors. Zh. Eksp. Teor Fiz. 34, 735-739 (1958) 
[Sov Phys. JETP 7, 505-508 (1958)]. 

[35] Eilenberger, G. Transformation of Gor'kov 's equation for type 2 superconductors into transport-like 
equations. Z. Phys. 214, 195-213 (1968). 

[36] Usadel, K.D. Generalized diffusion equation for superconducting alloys. Phys. Rev. Lett. 25, 507-509 

10 



(1970). 

[37] Bergeret, F.S. & Volkov, A.F. & Efetov, K.B. Josephson current in superconductor-ferromagnet struc- 
tures with a nonhomogenous magnetization. Phys. Rev. B 64, 134506 (2001). 

[38] Zaitsev, A.V. Quasiclassical equations of the theory of superconductivity for contiguous metals and 
the properties of constricted microcontacts. Zh. Eksp. Tear. Fiz. 86, 1742 (1984) [Sov. Phys. JETP 59, 
1015-1024 (1984)] 

[39] Kupriyanov, M.Yu. & Lukichev, V.F. Influence of boundary transparency on the critical current of 
"dirty" SS'S structures. Zh. Eksp. Tear. Fiz. 94, 139 (1988) {Sov. Phys. JETP 61, 1163-1168 (1988)] 

[40] Blonder, G.E. Tinkham, M. Klapwijk, T.M. Transition from metallic to tunneUng regimes in super- 
conducting microconstrictions: Excess current, charge imbalance, and supercurrent conversion. Phys. 
Rev. B 25, 4515(1982). 



11 



5 10 15 20 25 30 1 2 3 4 5 6 

Magnetic field (T) IVlagnetic field (T) 

c I I d — ■ — 1 — ■ — I — ■ — I — ■ — I — ■ — 1 — ■ — I — • — I — ■ — I — ■ — I 




Temperature (K) Landau level index (n) 



FIG. 2: Magnetoresistance oscillations of the Bi2Te3 surface states, (a) Angle dependence of the 
Shubnikov-de Haas oscillations, at 4.2 K, (top to bottom: the angle between the magnetic field and the sur- 
face normal = 0° to 60° in steps of 10°) after linear background subtraction. The square markers represent 
the expected shift of the maxima in magnetic field for a two dimensional system, given by B± = B cos(</>), 
indicating that the oscillations arise from surface states, (b) Oscillation amplitude dependence on the mag- 
netic field (minima as black squares and maxima as red circles); linear backgrounds have been subtracted 
leaving the oscillating part. We deduce a Dingle temperature = 1.65 K from fitting the increase in 
oscillation amplitude with increasing field. Only the magnetoresistance in small fields is considered, where 
the peak splitting is negligible, (c) Oscillation amplitude dependence on temperature; linear backgrounds 
have been subtracted. The effective mass m* = 0.16 mo is estimated by fitting the decrease of the mag- 
netoresistance oscillations as function of temperature (peaks at least up to 9 T in both the conductance and 
resistance result in the same effective mass), (d) The 1/B values versus the Landau level index n intersect 
at n = —0.5, consistent with a half filled lowest Landau level as expected for a Dirac cone. 
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FIG. 3: Josephson effects, (a) Current voltage characteristics of a Nb-Bi2Te3-Nb Josephson junction 
(length 1=50 nm) upon increased 10.0 Ghz RF irradiation power (curves are shifted for clarity). The tem- 
perature is 1.6 K. Clear Shapiro steps are observed at multiples ofV = ^/rf = 20.7 jiV. (b) Differential 
resistance, dV/ dl, plotted on a color scale as a function of the bias current, /, and the microwave excitation 
power, Prf- The numbers correspond to the n-th order Shapiro step, (c) Power dependence of the first three 
Shapiro steps. Shapiro steps are resolved in the junctions with length up to 150 nm. (d) Critical current 



sine ( ^ 



dependence on the magnetic field fitted with Ic = Iq 
flux and $o the flux quantum. The temperature is 260 mK. 



, /q the critical current at zero field, $ the 



13 



< 



10 



:-E 1 
6 



- Dirty limit, Usadel 
-Clean limit, Eilenberger 
Measurements 




50 100 150 200 250 300 
Junction length (nm) 



< 



o 



O 



Dirty limit, Usadel 
Clean limit, Eilenberger 
Measurements 




1 2 3 4 5 

Temperature(K) 



FIG. 4: Temperature and length dependence of the critical current; demonstration of the ballistic 
nature of the junctions, (a) Scaling of the junction critical current with the electrode separation length at 
1.6 K. (b) Temperature dependence of the critical current of the 50 nm junction. The measured data (red 
squares) can only be consistently fitted in both cases with Eilenberger theory for ballistic junctions (dashed 
black curves), while Usadel theory for diffusive junctions (dotted blue curves) cannot provide a good fit 
to the temperature dependence of the critical current. See the Supplementary Information for details. The 
obtained clean limit coherence length ^ = 75 nm at 1.6 K. 
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